Fragmentation in Peripheral Heavy-Ion Collisions: from Neck Emission to 

Spectator Decays 



o 
o 

(N 

c 
a 
>—> 

oo 

(N 



> 
OO 



o 
m 
o 

■ft 

=5 



X 



J. Lukasik, 1 - 10 G. Auger, 2 M.L. Begemann-Blaich, 1 N. Bellaize, 4 R. Bittiger, 1 F. Bocage, 4 
B. Borderie, 3 R. Bougault, 4 B. Bouriquet, 2 J.L. Charvet, 5 A. Chbihi, 2 R. Dayras, 5 D. Durand 4 
J.D. Frankland, 2 E. Galichet 3 ' 11 D. Gourio, 1 D. Guinet, 6 S. Hudan, 2 B. Hurst, 4 P. Lautesse, 6 F. Lavaud, 3 
A. Le Fevre, 1 R. Legrain, 5 '0 0. Lopez, 4 U. Lynen, 1 W.F.J. Miiller, 1 L. Nalpas, 5 H. Orth, 1 E. Plagnol, 3 
E. Rosato, 7 A. Saija, 8 C. Schwarz, 1 C. Sfienti, 1 J.C. Steckmeyer, 4 B. Tamain, 4 W. Trautmann, 1 
A. Trzcinski, 9 K. Turzo, 1 E. Vient, 4 M. Vigilante, 7 C. Volant, 5 B. Zwieglinski, 9 and A.S. Botvina 1 ' 12 

(The INDRA and ALADIN Collaborations) 

I Gesellschaft fur Schwerionenforschung mbH, D-64291 Darmstadt, Germany 

2 G ANIL, CEA et IN2P3-CNRS, F-14076 Caen, France 
5 1 Institut de Physique Nucleaire, IN2P3-CNRS et Universite, F-91406 Orsay, France 
4 LPC, IN2P3-CNRS, ISMRA et Universite, F-14050 Caen, France 
5 DAPNIA/SPhN, CEA/Saclay, F-91191 Gif sur Yvette, France 
6 Institut de Physique Nucleaire, IN2P3-CNRS et Universite, F-69622 Villeurbanne, France 
7 Dipartimento di Scienze Fisiche e Sezione INFN, Univ. Federico II, 1-80126 Napoli, Italy 
8 Dipartimento di Fisica dell' Universita and INFN, 1-95129 Catania, Italy 
9 A. Soltan Institute for Nuclear Studies, Pl-00681 Warsaw, Poland 
10 H. Niewodniczanski Institute of Nuclear Physics, Pl-31342 Krakow, Poland 

II Conservatoire National des Arts et Metiers, F75141 Paris cedex 03, France 

12 Institute for Nuclear Research, 117312 Moscow, Russia 
(Dated: February 8, 2008) 

Invariant cross sections of intermediate mass fragments in peripheral collisions of 197 Au on 197 Au 
at incident energies between 40 and 150 MeV per nucleon have been measured with the 4-7T multi- 
detector INDRA. The maximum of the fragment production is located near mid-rapidity at the 
lower energies and moves gradually towards the projectile and target rapidities as the energy is 
increased. Schematic calculations within an extended Goldhaber model suggest that the observed 
cross-section distributions and their evolution with energy are predominantly the result of the clus- 
tering requirement for the emerging fragments and of their Coulomb repulsion from the projectile 
and target residues. The quantitative comparison with transverse energy spectra and fragment 
charge distributions emphasizes the role of hard scattered nucleons in the fragmentation process. 
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The production of nuclear fragments is a decay mode 
that persists up to the highest bombarding energies in 
heavy-ion collisions. Its characteristics, however, are 
considerably different in the Fermi energy regime and 
in the relativistic regime. At bombarding energies near 
the Fermi energy, the fragment distributions are concen- 
trated at mid-rapidity, neck breaking mechanisms have 
been suggested by various transport simulations, and 
strong Coulomb effects are generated by the presence of 
slowly moving heavy residues [JSHSHSESHES 

ElElElEi. 

At relativistic energies, the fragment production in pe- 
ripheral collisions is concentrated at projectile and target 
rapidities El El [13 

with invariance properties that sug- 
gest a high degree of equilibration for the decaying spec- 
tator systems 0, 0, |2(j ■ Because of their prevailing 
statistical character these fragmentations are often seen 
distinguished from the dynamical processes at interme- 
diate energies. 
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The present study is directed to the transition region 
between these seemingly so different regimes with the 
aim to identify the basic ingredients that govern the evo- 
lution of the fragment production mechanisms. As the 
main result, we find that the same criteria are at work 
throughout the covered energy range and that, of these, 
the clustering (coalescence) criterion is essential for pro- 
ducing the observed features. 

The experiments were performed at the GSI Darm- 
stadt with beams of 197 Au, delivered by the heavy-ion 
synchrotron SIS and directed onto 197 Au targets of 2- 
mg/cm 2 areal thickness placed inside the INDRA mul- 
tidetector |2lj . The incident energies varied between 40 
to 150 MeV per nucleon, a range of energies that cor- 
responds to relative velocities between once and twice 
the Fermi value. Further details of the experimental and 
calibration procedures may be found in Refs. [22ll23|| . 

As described previously, the total transverse energy 
E] 2 of light charged particles [Z < 2) has been used as 
an impact parameter selector [U The measured E] 2 
spectra scale linearly with the incident energy, and the 
scaled spectra coincide. The relation between E 12 and 
the reduced impact parameter 6/6 max has been obtained 
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by applying the geometrical prescription of Ref . [2J| . It 
is linear in very good approximation, where &/6 m ax de- 
creases with increasing E] 2 . Eight finite impact param- 
eter bins were generated, with the most central bin 8 
covering impact parameters up to 5% of & max - The re- 
maining part of the E 12 spectrum has been divided into 
7 bins of equal width, corresponding to 7 bins of approx- 
imately equal width in &, of which bin 1 contains the 
most peripheral events. An event was registered if at 
least three detectors had recorded a hit. 
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FIG. 1: Invariant cross section distributions as a function 
of transverse velocity y/3± and rapidity y for fragments with 
Z — 3, 6, 9 from the most peripheral collisions (bin 1) of 
197 Au + 197 Au at E/A = 40, 60, 80, 100 and 150 MeV, as 
indicated. The cross sections are normalized relative to each 
other, separately for each fragment species. Near the target 
rapidity they are affected by the thresholds for fragment iden- 
tification. The arrows indicate the target, center-of-mass and 
projectile rapidities (from left to right). 

Invariant cross section distributions for lithium, car- 
bon, and fluorine fragments emitted in peripheral 197 Au 
+ 197 Au collisions are shown in Fig. ^ They illustrate 
the characteristic evolution of the emission pattern as 
the incident energy is raised from 40 up to 150 MeV per 
nucleoli. 

The Coulomb repulsion from the projectile and target 
residues that survive the collision at peripheral impact 
parameters is evident in all distributions. The distri- 
bution of intermediate-mass fragments, however, is not 
isotropic with respect to these residues. They are emitted 
with a preference toward mid-rapidity which produces 
a pronounced forward-backward asymmetry relative to 
the rapidities of the projectile and target residues. The 
resulting cross-section maximum at mid-rapidity at the 



lower energies gives way to a saddle-type structure at 
mid-rapidity at the higher energies. Here the fragment 
production is increasingly concentrated at the projectile 
and target rapidities, a trend that is more pronounced 
at even higher incident energies and generally character- 
istic for the relativistic regime [3. At a given energy, 
the probability for forming fragments at mid-rapidity de- 
creases rapidly with increasing Z . 

The kinetic properties of the mid-rapidity fragments in 
peripheral collisions, in spite of this apparent evolution, 
have been found rather constant throughout the covered 
range of projectile energies j2^|. The transverse velocities 
are invariant with respect to both, the incident energy 
and the fragment Z . This property and the large value 
(E t ) ~ 30 MeV of the corresponding mean transverse 
energies were suggested to be the result of the combined 
action of the intrinsic Fermi motion, the acquired trans- 
verse energies of scattered nucleons and of the Coulomb 
repulsion from the heavy residues. 

To obtain more quantitative predictions, a Monte- 
Carlo procedure based on these considerations has been 
developed. The original idea of Goldhaber [25| of tak- 
ing the fragment momentum as the sum of the momenta 
of the constituent nucleons picked randomly from a sin- 
gle Fermi sphere has been extended by including also 
the Pauli-allowed distribution of hard scattered nucle- 
ons. The corresponding calculations proceed in three 
main steps: 

(i) The momenta of the A nucleons of a fragment 
with mass number A are randomly picked from 2 Fermi 
spheres (pp = 265 MeV/c), separated by a relative mo- 
mentum per nucleon calculated at the distance of closest 
approach of the two nuclei colliding at a given impact pa- 
rameter, and also from the momentum distribution gen- 
erated in hard two-body collisions of nucleons from the 
projectile and target. Isotropic scattering has been as- 
sumed for these collisions, and they were requested to 
obey the Pauli principle. The Pauli-allowed momentum 
space changes with increasing incident energy and widens 
both toward lower and higher transverse momenta. In 
the calculations, a fixed impact parameter of 10 fm was 
used, and the number of hard scattered nucleons to be in- 
corporated in a fragment was taken as a parameter. The 
results presented here were obtained assuming a binomial 
distribution of the number of hard scattered nucleons in a 
fragment with mean values that were found to be rather 
small (see below). 

(ii) A clustering criterion has been implemented by 
randomly accepting or rejecting the trial clusters accord- 
ing to the probability distribution of the mean squared 
relative momentum, ^2 i: £j(pi—Pj) 2 / (A(A—1)), with that 
of A nucleons drawn randomly from a single Fermi sphere 
of the radius of pp = 265 MeV/c. Examples of these ref- 
erence distributions are shown in the upper part of Fig. [21 
for fragments containing 2, 7, 14, and 28 nucleons. They 
are characterized by a mean value (p^ei) anc ^ a variance 
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FIG. 2: Top panel: Probability distributions of the mean 
squared relative momentum of nucleon pairs picked randomly 
from a single Fermi sphere for fragments with A = 2 (solid 
line) and Gaussian approximations of the distributions for A 
= 7, 14 and 28 (dashed, dotted and dash-dotted lines, respec- 
tively). Bottom panel: Initial spatial conditions for Coulomb 
trajectory calculations in the reaction plane, x vs z. The lines 
represent the trajectories of the projectile and the target at 
incident energies E/A — 40 to 150 MeV and for a fixed impact 
parameter of 10 fm. The dots placed on the trajectories mark 
the positions of the centers of the nuclei at the distance of 
closest approach (gray) and in the exit-channel sticking con- 
figuration (black). The arrows indicate the directions of the 
ion motion. 
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The mean value of the distributions is constant while 
the variances decrease with increasing fragment mass A. 
The asymmetry apparent for lighter fragments disap- 
pears for heavier ones and, for A > 7, the distributions 
were approximated by Gaussians. The clustering crite- 
rion strongly suppresses the formation of fragments with 
large relative momenta of constituent nucleons. 

The fragments selected by random picking of the mo- 
menta from a single Fermi sphere are not entirely cold 
(see e.g. 26]), and neither are the fragments defined 
with the above criterion. Temperatures deduced by fit- 
ting a Fermi distribution to the distribution of nucleon 



momenta in the center of mass of a fragment drops, ap- 
proximately like 1/VA, from about 5 MeV for A = 10 to 
near for very large clusters. The sequential decay of ex- 
cited clusters is expected to modify the fragment yields 
and kinetic energies, but not to an extent that would 
substantially modify the present conclusions. 
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FIG. 3: Invariant cross sections of 7 Li fragments as a function 
of transverse velocity 7/3x and rapidity y, as obtained from 
the model calculations for peripheral 197 Au + 197 Au collisions 
at E/A = 40 MeV (top part) and 150 MeV (bottom). The 
left and the right columns show the initial and asymptotic 
distributions, i.e. before and after the Coulomb trajectory 
calculations, respectively. The results without nucleon mix- 
ing, without the clustering criterion, and of the full calculation 
are shown in the top, middle, and bottom rows for both en- 
ergies, respectively. The arrows have the same meaning as in 
Fig- IH 

The specific clustering criterion is adapted from the 
original Goldhaber formulation of the fragmentation 
problem by generalizing it from a single Fermi sphere 
to the three distributions considered here. It differs from 
standard coalescence criteria based on a fixed momentum 
radius and reduces the number of free parameters of the 
model as it depends only on pp. 

(iii) Accepted fragments have been placed, in coordi- 
nate space, in between the two residues in an aligned 
sticking configuration which served as the starting point 
for 3-body Coulomb trajectory calculations (Fig. |5J bot- 
tom panel). The initial momentum of the produced 
fragment is determined by steps (i) and (ii), while the 
starting momenta of the heavy fragments were generated 
with two-body trajectory calculations for the entrance 
channel, using Coulomb and proximity forces |2?|. The 
residue trajectories remain virtually undisturbed in the 
exit channel. This simplified reaction picture assumes 
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that the fragment is formed in the contact zone of the 
colliding nuclei and released as soon as the remaining 
residues separate by more than the diameter of the frag- 
ment which is assumed to be spherical. The effect of 
additional charged particles is ignored since their multi- 
plicity in peripheral collisions is small. 

The results obtained with this procedure for 40 and 150 
MeV per nucleon, the lowest and highest bombarding en- 
ergy included in this study, are shown in Fig. |3| for the 
case of fragments with mass number A = 7. Three dif- 
ferent conditions were chosen for the calculations which 
all have the assumption in common that the number of 
hard scattered nucleons in a fragment is given by a bino- 
mial distribution with the mean value of 1. The resulting 
cross section distributions are rather similar at the lowest 
energy but become distinctively different as the energy is 
raised. 
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FIG. 4: Experimental (symbols) and calculated (lines) 
transverse-energy spectra of 7 Li (left panel), 13 C (middle), 
and 19 F fragments (right) emitted at mid-rapidity in periph- 
eral 197 Au + 19r Au collisions at five incident energies, as in- 
dicated. Solid, dashed and dotted lines correspond to the cal- 
culations with on average one, zero or two scattered nucleons 
in a fragment, respectively. In each panel, the experimental 
spectra are displaced vertically by consecutive factors of 30, 
and the calculated spectra are individually normalized rela- 
tive to the corresponding measurement. 

If all the non-scattered nucleons forming a fragment 
are exclusively taken from only one of the collision part- 
ners, the mid-rapidity region cannot be populated at the 
higher bombarding energies (top row of Fig. [3] for 150 
AMeV). At the lower energies (top row for 40 AMeV) 
these evaporation-like events are copiously found at mid- 
rapidity because the Coulomb maxima, forming the char- 
acteristic rings in the two-dimensional representation, 
intersect. If the clustering criterion is dropped, (mid- 
dle rows for both energies), the fragment distributions 
at high incident energies populate nearly homogeneously 
the full range from projectile to target rapidities, again 
in contrast to the observation (cf. Fig. [IJ. Both ingre- 
dients, the mixing of target and projectile nucleons (see 
also |6|) and the clustering criterion, are indispensable if 



more realistic distributions are to be generated that re- 
semble those observed (bottom rows for both energies). 
The Coulomb field and its interplay with the initial mo- 
menta has a strong effect of rearranging and focusing 
the asymptotic fragments (cf. Ref. [111). Their distribu- 
tion extends to the projectile and target rapidities even 
though all fragments are initially placed in the aligned 
sticking configuration in between the two residues (see 
(iii) above). These specific initial conditions naturally 
produce the forward-backward asymmetry with respect 
to the heavy-residue rapidities. 

For a more quantitative comparison, measured and cal- 
culated transverse-energy spectra for mid-rapidity frag- 
ments with Z = 3, 6, 9 are shown in Fig. 21 Bins cen- 
tered at mid-rapidity with widths of 25% of the projectile 
rapidity y p were chosen for the selection. The exponen- 
tial spectra shapes correspond to normal statistical emis- 
sions. Only the high value of the mean transverse energy 
(« 30 MeV) has been found difficult to be explained 
within thermal models, as discussed in Ref. |22|. The 
shapes of these spectra, including their invariance with 
bombarding energy and Z, are rather well reproduced if 
on average one scattered nucleon from the Pauli-allowed 
distribution of hard scattered nucleons is included in the 
fragment (solid lines) . If scattered nucleons are excluded 
(dashed lines) or if their mean number is increased to 2 
(dotted lines), the resulting spectra deviate significantly 
from the experiment, in particular for the lighter frag- 
ments. 

A small number of scattered nucleons in a fragment 
is consistent with the small number of primary nucleon- 
nucleon collisions expected for the very peripheral col- 
lisions selected with bin 1. The increasingly harder 
transverse-energy spectra in more central collisions |22l | 
require more scattered nucleons to be included, in par- 
ticular at the higher incident energies. For 100 and 150 
MeV per nucleon, the lithium spectra of bins 2 and 3, 
e.g., are well described if their mean number in a frag- 
ment is increased to 1.5 and 2, respectively. 

The general trends of the experimental cross section 
distributions as a function of the fragment Z are also 
rather well reproduced. This is demonstrated in Fig. [5] 
for the range of mid-rapidity fragments with 3<Z<10. 
The Z distributions are approximately exponential func- 
tions with slopes that increase monotonically with the 
incident energy. The calculations have been performed 
for masses A = 2Z + 1 and for, on average, one scat- 
tered nucleon in a fragment. The Z dependence of the 
differential multiplicities is nearly perfectly reproduced 
for 100 MeV per nucleon but the variation of the slopes 
with incident energy is slightly overpredicted. With the 
chosen normalization at Z = 6, this has the effect that 
the yields of the heavier fragments are underestimated at 
the highest and overestimated at the lower bombarding 
energies. The yields of lithium fragments are also too 
low by about a factor of two at all energies except 150 
MeV per nucleon, where the measured yields are lower. 
Overall, however, the agreement is very satisfactory. 
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FIG. 5: Experimental (symbols) and calculated (lines) mean 
multiplicity of mid-rapidity fragments in peripheral 197 Au + 
197 Au collisions (bin 1) as a function of Z. The calculated 
spectra are normalized at Z — 6, separately for each of the 
five bombarding energies E/A = 40, 60, 80, 100, 150 MeV 
(from top to bottom). 



From the good qualitative (Figs. n an d|HJ and quan- 
titative agreement (Figs. 0] and EJ) we may thus con- 
clude that the presented extended Goldhaber scenario 
contains the essential ingredients of the fragmentation 
mechanism in peripheral collisions. Of these, the clus- 
tering criterion is responsible for the evolution of the 
cross section distributions from mid-rapidity to the spec- 



tator rapidities with increasing projectile energy. It is 
furthermore important to permit the mixing of projectile 
and target nucleons and to include hard-scattered nucle- 
ons from Pauli-allowed collisions into the produced frag- 
ments. This latter observation even suggests the hypoth- 
esis that the hard-scattered nucleons may be required 
for initiating the fragment production at mid-rapidity in 
peripheral collisions. The correct prediction of the role 
of scattered nucleons may therefore be considered as an 
important test for any dynamical model to describe the 
mid-rapidity emission. It is, finally, remarkable how far 
the fragment distributions extend into the projectile and 
target rapidity regions, as a consequence of the combined 
action of the initial momenta and of the Coulomb forces 
giving rise to the characteristic ring structures. Appar- 
ently, fragments observed there do not necessarily have 
to be considered as emitted from the excited residues but 
may to a large part originate from the contact zone, or 
neck region, formed during the reaction. 
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